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ABSTRACT 

Context. Our comprehension of stellar evolution on the AGB still faces many difficulties. To improve on this, a quantified understand- 
ing of large- amplitude pulsator atmospheres and interpretation in terms of their fundamental stellar parameters are essential. 
Aims. We wish to evaluate the effectiveness of the recently released CODEX dynamical model atmospheres in representing M-type 
Mira variables through a confrontation with the time-resolved spectro-photometric and interferometric PTI data set of TU And. 
Methods. We calibrated the interferometric K-band time series to high precision. This results in 50 nights of observations, covering 
8 subsequent pulsation cycles. At each phase, the flux at 2.2 yum is obtained, along with the spectral shape and visibility points in 5 
channels across the K-band. We compared the data set to the relevant dynamical, self-excited CODEX models. 
Results. Both spectrum and visibilities are consistently reproduced at visual minimum phases. Near maximum, our observations show 
that the current models predict a photosphere that is too compact and hot, and we find that the extended atmosphere lacks HiO opacity. 
Since coverage in model parameter space is currently poor, more models are needed to make firm conclusions on the cause of the 
discrepancies. We argue that for TU And, the discrepancy could be lifted by adopting a lower value of the mixing length parameter 
combined with an increase in the stellar mass and/or a decrease in metallicity, but this requires the release of an extended model grid. 
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1. Introduction 



In sp ite of their great astrophysical importance dMuzzin et alJ 
2009), many aspects of asym ptotic giant branch ( AGB) stars re- 
main poorly understood (e.g. iKerschbaum et alj|201 U and ref- 
erences therein). Among the AGB stars, the luminous Mira 
variables might be the most enigmatic with their long-period, 
large-amplitude pulsations and high mass-loss rates. In particu- 
lar our understanding of their atmospheric structure, which is the 
link between enriched stellar interior and evolution-dominating 
wind, is at a turning point, thanks to advances on both the obser- 
vational side (IR interferometry), and the modelling side with a 
new generation of self-excited atmosphere models. 

The public rel ease of four CODE X model series for M-type 
Mira variables bv llreland et alJ d201ll hereafter ISW1 1), greatly 
increases the opportunities for the general community to assess 
the performance of such dynamic models. Here, we present the 
first confrontation of the CODEX models (see Sect. O to a fully 
phase-resolved 2 yum spectro-photometric and interferometric 
dataset covering eight pulsation cycles, as obtained with the 
Palomar Testbed Interferometer on the Mira TU And. This star 
has a spectral type M5IIIe, a pro genitor mass M - 1 . 15 - 1 .4 M , 
metallicity Z = 0.004 - 0.02 dMennessier et alJ l200lh . and a 
GCVS period P-317 d (ISamus et al.ll2009h . 

The most intriguing result of this comparison is that there 
seems to be a systematic discrepancy between the models and 
our observations around visual maximum: the star appears too 
hot, and its outer atmosphere too devoid of water. An identifica- 
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tion of the source of this discrepancy is hampered by the limited 
coverage in parameter space of the models, but we discuss pos- 
sible avenues of further research, and stress the need for larger 
published model grids and for observational campaigns covering 
the visual maximum phase in detail. 

2. Observations 

TU And was extensively observed as part of the large PTI Mira 
programme (iThompsonl |2002|) during 60 nights from 1999 to 
2005, over eight pulsation cycles. The publicly available PTI 
data archive provides both squared visibility amplitudes (here- 
after visibilities) and integrated photon counts. We describe the 
data reduction and calibration of these observables below. 



2.1. Interferometry 

The now dismounted Palomar Testbed Interferometer (PTI) 
was a Y-shaped, long-baseline near-IR interferometer located 
at Palomar Observatory and operated by JPL/Caltech. It con- 
sisted of three 40 cm apertures with separations of 85-1 10m 
and used in pairwise combination (see IColavita et alJIl 999L for 
a detailed description). PTI thus measured only one interfero- 
metric observable per spectral channel: the square of the "fringe 
contrast" or visibility amplitude. Visibilities were recorded in 
five channels within the K-band (2.0-2.4 /im). Processing of the 
raw data (correcting for detector biases and read-out noise, and 
averaging the 125s scans into fringe amplitudes) was automati- 
cally done on site, following procedures described in IColavita! 
(1 19991) . Calibrating for the system visibility, i.e. the response 
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of the measurement system (interferometer + atmosphere) to 
a poi nt source, was done using the nbCalib tool dBoden et alJ 
119981) provided by the NExScI which uses a weighted-mean 
formalism to estimate the system visibility at the time and po- 
sition of the science observation from all calibrator measure- 
ment s that fit the user-def ined temporal and spatial constraints 
(see Ivan Belle et alii 1999L for more details). To ensure calibra- 
tion stability and consistency, essential for making full use of 
this data set as an interferometric time series, several calibrators 
were observed each night and were intercompared by nbCalib. 

The stellar parameters of our calibrators w ere determined by 
fitting a reddened Kurucz model atmosphere dKuruczll 19931) to 
high-quality archival photometry. The errors were estimated us- 
ing Monte Carlo simulations. Table Q] (only available electroni- 
cally) lists the calibrators and their adopted parameters. All cal- 
ibrators have angular diameters well below the PTI resolution 
limit of 1 mas (Ivan Belle & van BeUl]2005) and are, on average, 
10° away from the science target. 

The calibrated visibilities of TU And from 50 nights were 
retained, having a median relative precision ranging from 5% at 
2.2 ytzm, 7.5% at 2.4 ^m and 10% at 2.0 yum. 

2.2. Spectro-photometry 

Together with the uncalibrated visibilities, PTI delivered the raw 
photon counts of each observation, typically in blocks of five ob- 
servations taken within minutes and interleaving science targets 
with calibrators. We derived the spectral shape over the K-band 
and flux-calibrated the 2.2 /mi channel with a bootstrap method. 

The spectral shape calibration procedure consists of divid- 
ing the science-target normalised photon counts by the ones 
of the calibrator and multiplying the result with the intrinsic 
spectral shape of the calibrator. To obtain an absolute flux in 
the 2.2 fim channel, one calibrator measurement on either side 
of the science observation with respect to airmass is required 
to solve for the atmospheric extinction per unit airmass A v = 
(log[/2/^i] + log[5 \/S2])/{mi - fn{) and the instrument spectral 
response R v = S x exp(— A v mi)), with the intrinsic and 5, 
the measured flux, and m, the airmass, of the z'th calibrator. The 
calibrated flux is then calculated as 7 JC , = S Jf , ■//?,. x exp(A v m JC! ), 
where S denotes the measured photon counts and m JC , the air- 
mass of the science target. For each observing block, the above 
procedure is applied 5000 times to a randomly picked science 
observation, corresponding random calibrator observations and 
calibrator intrinsic fluxes. Finally, the distributions of all observ- 
ing blocks within a single night are co-added, and the average 
and standard deviation are taken from the resulting normal dis- 
tribution, as the final calibrated flux and its error. 

To obtain the calibrator intrinsic fluxes, the PTI optical con- 
figuration was taken into account: after beam combination and 
dispersion by a prism, the wavefront was s ent through a broad - 
band K Barr filter onto the detector pixels (IColavita et al. 1999). 
Each of the Kurucz models in the Monte Carlo distributions ob- 
tained with the spectral energy distribution fitting, is converted 
into the effective flux I in each of the PTI channels as: 
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with M the Kurucz model, F the filter profile, G a Gaussian of 
FWHM = 0.7 fim (M. Colavita, private communication), and A/, 
A u the edges of the spectral channel under consideration. 

Finally, we constructed a visual light curve by retrieving ob- 
served V-band photometry from the AAVSO and AFOEV. Visual 



phases (hereafter simply phases) were computed by fitting a sine 
to ten cycles around the times of the PTI observations. 

3. The CODEX models 



The CODEX atmosphere models dlreland et all 120081) are the 
most advanced publicly available dynamical models for M-type 
Mira variables. The pulsations in these models are self-excited 
and limited in amplitude by the non-linear effects of shocks in 
the outer atmosphere and by turbulent viscosity. For each model 
series, a few typical cycles were selected for detailed radiative 
transfer calculations with an opacity-sampling method in LTE 
with a phase-resolution of ~ 0.1. This results in the centre-to- 
limb variation (CLV) and spectrum, with AA = 0.0002 fim be- 
tween A = 0.45 - 2.50/mi. The current model series are made 
to match o Cet, R Leo, and R Cas, although even these well- 
observed Miras have large uncertainties on their basic parame- 
ters (e.g. the 2 very different models for R Cas). They were pre- 
sented to the community and compared to photometric and in- 
terferometric observations in ISW1 1. The model parameters are 
mass M, luminosity L, metallicity Z, mixing length parameter 
a m , and turbulent viscosity parameter a v . For a detailed justifi- 
cation of the current set of parameter values, we refer to ISW1 1 . 

With a slightly meandering (+ 5 d) pulsation period of 
-317 d, TU And is compared to the R52 (M = 1 . 1 M , L = 
5200 L Q , Z = 0.02, a m = 3.5, a v = 0.25) and o54 (M - 1.1 M Q , 
L = 5400 L , Z = 0.02, a m = 3.5, a v = 0.25) model series, 
which have derived periods of 307 and 330 d, respectively. 

When comparing a model series to observations, the only 
free parameter is the distance, or equivalently, the phase- 
averaged Rosseland angular diameter Or. Considering the angu- 
lar resolution of PTI, we converted the model intensities to syn- 
thetic observables for Or e [1.5,4.5] mas, in steps of 0.05 mas. 
Converting the model intensities to the spectro-photometric ob- 
servables was d one with Eq.[D The squ ared visibilities are com- 
puted following lWittkowski etaTI d2004l) : 
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with F, G, Ai, and A u as in Eq. [1] and 
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the monochromatic Hankel transform of the intensity profile I A 
at the baseline B and wavelength A. 

4. Results 

Figure [1] shows the resulting photometric time series of TU 
And as a function of phase, with two consecutive cycles of the 
R52 model series repeatedly overplotted. The two other avail- 
able cycles of this series give qualitatively identical results, so 
are not shown. The Rosseland angular diameter is chosen to 
match the absolute flux at 2.2 fim, with a resulting value of 
6 R = 2.20 + 0.15 mas or d = 1.13 + 0.07 kpc. Although the 
error on the absolute calibration in this channel is rather large, 
the contamination from upper atmospheric molecular layers at 
these wavelengths is smallest. This channel thus provides the 
best view at the continuum photosphere in the spectral region 
where the bulk of the luminosity is emitted. Moreover, fixing the 
angular diameter like this gives a consistent picture in both the 
spectro-photometry and the interferometry. 
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Fig. 1. Spectro-photometry. The two upper pan- 
els show the variation in the spectral shape with 
respect to the 2.2 fim channel, the third panel 
shows the absolute flux in the latter, and the 
lower panel the visual flux in magnitudes. In 
each panel two cycles of the R52 model se- 
ries are consecutively overplotted as the solid 
line. On the left, the cycles are shown as a func- 
tion of visual phase while on the right they are 



The R52 model series shows reasonable agreement with the 
observed visual light curve in the lower panel, considering that 
(1) the amplitude was not tuned to match this light curve and 
2) non-LTE effects play an important role at these wavelengths 
Ireland et al . 2008), causing an increase in visual brightness of 
0.5 to 1.0 magnitude. In case the higher value of this correction 
is valid, the perceived difference between model and observa- 
tions could point to the presence of an excess in visual light at 
maximum phases together with a too large model amplitude. 

The two upper panels show the variation in the spectral shape 
with time. A clear discrepancy occurs around times of visual 
maximum, then the model predicts too much flux. Moreover, 
the amplitude of this excess shows an interesting wavelength- 
dependence: the strongest effect is observed at 2.0 yum, a smaller 
discrepancy at 2. 1 and 2.4 jjm and a good agreement at 2.3 //m. 

The visibility time series is presented in Fig. [2] overplotted 
with the same cycles of the R52 model series. Each panel dis- 
plays the visibilities at a given phase, using different symbols 
for observations coming from different cycles (at least two per 
panel). The model reproduces the observations well at phases 
0.3-0.6, but not around visual maximum. Comparison with other 
cycles does not give better results. In fact, this reveals that 
the model has larger cycle-to-cycle variations at pre-maximum 
phases than are seen in the observations. A closer look at the 
wavelength dependence shows that the discrepancy is least pro- 
nounced at 2.3 and 2.4 yum. Between 2.0 and 2.2 yum, the model 
shows very little wavelength-dependence, unlike the observa- 
tions which show a clear decrease in visibility towards 2.0 fim. 

This observation is supported by the more quantitative as- 
sessment presented in Figs. [3] and [4] (available electronically 
only). Epochs with observations on at least two different base- 
lines were selected, to resolve the degeneracy between CLV ba- 
sic size and shape. Each of these data sets with good uv-coverage 
was compared to each model, independently of model phase and 
for a range of 6r. The figures show the resulting^- 2 -maps. In all 
panels a curved valley shows up: every data set can be fitted with 
a model of any phase, albeit with a significantly larger diameter 
for visual maximum models. However, within these valleys a 
better x 2 is typically found with a model at the correct (photo- 
metric) diameter but with a phase nearer to visual minimum than 
the actual value of the observations. In conclusion, we find that 
the models around visual maximum should resemble those of vi- 
sual minimum more in CLV shape and wavelength-dependence. 

Finally, we also compared the observations to the similar o54 
model series, except for its luminosity. The results for o54 are 
alike and thus not shown here. The visibilities seem to match 
slightly better, but the light curve discrepancies become larger. 
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Fig. 2. Visibility time series. Each panel shows all squared visi- 
bilities as a function of spatial frequency, obtained at the visual 
phase given in the lower left corner of the panel. Different sym- 
bols are used for observations from different cycles. The same 
two model cycles as in Fig. Q] are overplotted as the full and 
dashed lines in each panel. 

5. Discussion 

llreland et ail (120081) demonstrate that the CODEX models have 
reached a level of maturity, concerning computational issues and 
input physics, that allows quantitative comparisons with obser- 
vations. Pulsations are self-excited and molecular layers are pro- 
duced naturally and in a physically self-consistent way. 

Nevertheless, our observations of TU And show a clear dis- 
crepancy with the models around phases of visual maximum, 
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Fig. 3. The x 2 maps of two data sets with a good uv-coverage 
(see text). The log y 2 is shown as a function of Rosseland an- 
gular diameter and phase, of the same two model cycles as in 
Figs. Q] and [2] The phase of the data set is given in the upper 
right corner of the panel and is plotted as the horizontal dotted 
line. The vertical dotted line marks the photometric diameter. 
The intersection of the dotted lines shows where the minimum 
X 1 should be if the model is correct. 



both in the spectrum and in the visibilities. The most natural in- 
terpretation of the wavelength-dependence of this discrepancy 
seems to be in terms of the shape of the H2O opacity curve, 
which shows an increase from 2.2 //m towards 2.0 and 2.4 fim. 
The formation of H2O layers depends sensitively on the density, 
temperature, and pressure stratification in the upper atmosphere. 
In the models, the conditions in the upper atmosphere are set 
by the dynamics. Shock fronts that emerge from the continuum 
photosphere between phases -0.3 and -0.1 and then propagate 
outwards lead to the formation of an H2O layer around phase 
0.2-0.3, which then strongly decreases in strength at phase 0.7- 
0.8. At (f>v = 0.5 both visibilities and spectra show good agree- 
ment. However, as demonstrated by Figs.[TJ|2] and [3] the near- 
maximum observations are represented better by model phases 
nearer to minimum, and thus with the presence of an H2O layer 
in the atmosphere. This interpretation is supported by the better 
agreement at 2.3 and 2.4 yum. There CO is an important source 
of opacity, in the photosphere and in the extended atmosphere, 
even at visual maximum. As a much more stable molecule, CO 
depends less on the exact conditions and partly hides the lack of 
an H2O molecular layer in the spectrum and visibilities. 

In retrospect, it is intere sting to note that the im age recon- 
struction of T Lep (P~380d, iLe Bouquin et al.ll2009h . made at a 
phase of 0.8-0.9, shows the clear signature of an extended molec- 
ular layer with a brightness ratio of ~10%, u nlike its closest 
model series o54, which supports our findings. Woodruff et al. 
(2009) found fair agreement between the wavelength depen- 
dence of their measured 1.1-3.8 /mi UDs and the predictions 
of the CODEX models, for several Miras at phases 0.5 to 0.8, 
but needed slightly larger model continuum diam eters than pre- 
dicted. On the other hand. lWittkowski et al.l d201 ll) have recently 
compared AMBER medium spectral resolution observations of 
four Miras to the CODEX models and found good agreement, 
although they compared the models to stars of a different period, 
hence luminosity. 

In summary, at visual maximum the current models predict 
a continuum photosphere that is too compact and hot as already 
suggested by ISW1 1, in combination with too small an H2O op- 
tical depth in the outer atmosphere, as found here. 

The basic question now is, can this be salvaged within the 
current modelling framework? Is it just a question of tuning the 
parameters to more realistic values? Following the discussion of 
ISW11, it seems difficult to find any combination of parame- 



ter values that fits all observables. Accepting the LMC period- 
luminosity relation and assuming that the effect of metallicity is 
insufficient (see ISW1 1), the only free parameters are M, a m , and 
a v . Since the last mainly influences the pulsation amplitude, but 
not the period or the 'parent star' radius, the above authors con- 
clude that the problem at hand is one of calibrating the mixing 
length parameter as a function of mass. Decreasing the rather 
high a m = 3.5 to a more 'realistic' value requires an increase 
in mass or a decrease in metallicity (albeit with a smaller de- 
pendency) to preserve the period, which would in case of o Cet 
be inconsistent with its kinematically determined age. TU And 
though, was classified as a member of a population with ini- 
tial masses in the 1.15 — 1.4 M range and Z e [0.004,0.02] 
dMennessier et alj |2001). which makes it a better candidate to 
test the above hypothesis. 

We can at this point not gauge the impact of using an a v 
which is not tuned to the amplitude of TU And. As a result, we 
feel that the model grid should span some scope in this parameter 
to allow more comparisons with 'small-amplitude' pulsators. 

Could either change in the model parameters also induce a 
stabilizing effect on the H2O layer, as required by our observa- 
tions, or is there a need to add more physics to the models? 

6. Conclusions 

The combination of a spectro-photometric and interferometric 
time series is a powerful tool for testing the validity of pulsat- 
ing AGB model atmospheres, even at low spectral resolution. 
PTI observed a large sample of M-type Miras and semi-regulars, 
several with a similar sampling to the one presented here, cov- 
ering periods between 150 and 470d. To fully exploit this rich 
database, better coverage in the model parameter space is ur- 
gently needed, even though it requires a large computational ef- 
fort. Only in this way can the model parameters be calibrated 
and (simply) linked to Mira observable properties. From the ob- 
servational point of view, interferometric studies with better uv- 
coverage and spectral resolution should focus on phases 0.0 to 
0.15, since none of the current model series predicts any ex- 
tended molecular atmosphere in that phase range. 
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Table 1. The derived calibrator parameters. 



Calibrator 


SpT 


7-eff (K) 


UD Diameter (mas) 


HD166 


KOV 


5540 ± 90 


0.61 


± 0.05 


HD1404 


A2V 


8840 ± 175 


0.47 


±0.02 


HD2628 


A7III 


7200 ± 50 


0.48 


±0.05 


HD3651 


KOV 


5250 ± 50 


0.74 


±0.06 


HD6920 


F8V 


6350 ± 200 


0.61 


±0.04 


HD7034 


FOV 


7600 ± 450 


0.49 


±0.04 


HD7229 


G9III 


5550 ± 150 


0.80 


± 0.12 


HD9712 


K1III 


4750 ± 240 


0.78 


± 0.20 




Fig. 4. Similar to Fig. [3] but for the four other data sets with good 
uv-coverage. The log x 1 is shown as a function of Rosseland 
angular diameter and phase, of the same two model cycles as in 
Figs. Q] and [2] The phase of the data set is given in the upper 
right corner of the panel and is plotted as the horizontal dotted 
line. The vertical dotted line marks the photometric diameter. 
The intersection of the dotted lines shows where the minimum 
X 1 should be if the model is correct. 



